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We address the concept of direct multiphoton multiple ionization in atoms exposed to intense, 
short wavelength radiation and explore the conditions under which such processes dominate over 
the sequential. Their contribution is shown to be quite robust, even under intensity fluctuations 
and interaction volume integration, and reasonable agreement with experimental data is also found. 



PACS numbers: 32.80.Rm, 32.80.Fb, 30.80.Hd, 42.50.Hz 
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The quest for direct multiple excitation/ionization of 
several electrons bound in atoms or molecules under the 
influence of intense laser radiation dates back to the early 
80 's Until the recent appearance of the FEL- based 

intense XUV and X-ray sources @, however, the hitherto 
available infrared and optical sources, although quite in- 
tense, proved inadequate for the task. Multiple ionization 
has of course been observed in a number of experiments 
but the mechanism has been the sequential strip- 
ping of electrons, beginning with the valence shell and 
proceeding inwards. 

The only notable exception to the sequential strip- 
ping, under long wavelength radiation (~ 800nm), came 
about with the advent of ultrashort (subpicosecond) TiSa 
sources, which led to the observation of non-sequential 
double ionization [?J . Still, the basic mechanism relies on 
a valence electron pulled out by the field, set into oscilla- 
tory motion thereby causing the ejection of a second, or 
perhaps third, electron by collision, as it returns towards 
the core; hence the term recollision Q. For this to be 
possible, the ponderomotive energy (cycle-averaged ki- 
netic energy) of a quasifree electron under the field must 
be larger than the binding energy of a second electron. 

The situation has now changed dramatically with the 
appearance of XUV to X-ray intense sources which has 
made feasible for the first time the observation of a num- 
ber of non-linear processes in that wavelength range 0- 
The decisive developments in that respect are the 
large peak intensity and sub-picosecond pulse duration. 
The latter is of central importance to our considerations 
in this work. An early, small scale so to speak, devel- 
opment in that direction has been the direct 2-photon 
double ionization of Helium which has grown into a sub- 
field with tens of theoretical [l2[ and a few experimental 
studies [HI , limited mainly by the present early stage of 
the sources. The chief difference between this process and 
its counterpart under long wavelength is that no recolli- 
sion is involved, as the ponderomotive energy is totally 
negligible. The two electrons are pulled out by the field 
while, electron-electron interaction although present, is 
not necessary; in contrast to single-photon double ion- 
ization. The information on this process, collected so 



far, can serve as a calibration for the larger scale gener- 
alization proposed in this paper. 

Although the idea, under suitable conditions, is appli- 
cable to essentially any atom, in the interest of providing 
a quantitative assessment of the underlying physics, we 
focus on the specific context of Neon under radiation of 
photon energy 93 eV (« 13.3 nm). In fact, some ex- 
perimental data have already appeared in the literature 
[141 ]. For this photon energy, even at a peak intensity 
10 18 W/cm 2 , the ponderomotive energy U p is about 10 
cV, which is much smaller than both the photon en- 
ergy and the binding energy of any electron in Neon, 
thus guaranteeing the validity of LOPT (Lowest non- 
vanishing Order of Perturbation Theory) 0, 0, [H| . In 
addition, for the notion of the cross section (generalized 
cross section of the appropriate order) to be valid, the 
pulse duration must be at least 10 cycles of the field. For 
photon frequencies of the order of 90-100 eV, even a pulse 
duration of lfs amply satisfies this condition. If the data 
are limited only to populations of the ionic species pro- 
duced in the process, a set of kinetic (rate) equations are 
sufficient for the interpretation, as well as for our pur- 
pose in this paper. This is a set of differential equations 
governing the evolution of the populations of the various 
ionic species during the pulse [H, [l5| . The complete set 
of such equations involving sequential, as well as direct 
processes from the neutral, leading to the ejection of up 
to 8 electrons are of the form 
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with 1 < j < 8 and <d[j] the discrete Heaviside function. 

The terms in the right hand side represent various pro- 
cesses contributing to the rate of change of the species 
whose time derivative appears on the left hand side. Thus 
Nj indicates the jth ionic species of charge while a 

(n) 

term like cr- kF n Nj represents an n-photon process lead- 
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ing from species j to species k, with the corresponding n- 
photon (generalized) cross section cr^. , where F(t) is the 
time-dependent photon flux in photons/cm 2 sec. These 
equations are to be solved under a pulse, as realistic as 
possible, dictated by the conditions of operation of the 
source; in this case the FEL. 

The simplest scenario consists of retaining only sequen- 
tial processes in the equations, which means that only the 
sequence of channels leading from ion jtoj+1, by succes- 
sive single electron ejections, is retained in the equations. 
At relatively low peak intensity and long pulse duration, 
these will be the dominant channels. As a point of cal- 
ibration for things to follow, we present in Fig. [1] the 
ion yields as a function of peak intensity, at the end of 
deterministic pulses of duration 30 fs [Fig. [TJa)] and 5 
fs [Fig. DIb)]. In both figures, the dashed lines repre- 
sent the yields for the sequential channels alone. These 
figures depict the typical single atom behavior, illustrat- 
ing the appearance and disappearance of ionic species as 
they give rise to higher ones with rising intensity. 
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FIG. 1: (Color online) Ionization of Ne at 93 eV under a 
deterministic pulse with duration (a) 30 fs and (b) 5 fs. Solu- 
tion of Eqs. fl]) in the presence of sequential channels alone 
(dashed lines) and with both direct and sequential processes 
included (solid lines). 



For Ne under 93 cV photons, however, an entirely new 
class of channels are energetically possible. These are 
direct, several electron multiphoton (of the appropriate 
order) processes, leading from the neutral to the corre- 
sponding ion. Specifically, 2 photons can eject 2 elec- 
trons leading to Ne +2 , 3 photons can eject 3 electrons 
leading to Ne +3 , etc., up to 6 photons leading directly 
to Ne +6 . These are higher ord er g eneralizations of 2- 
photon 2-electron ejection in He |12|. Moreover, we have 
included an 8-photon 7-electron transition Ne— > Ne +7 , 
and an 11-photon 8-electron transition Ne— > Ne +S . Di- 
rect n-photon m-electron ejection can in principle always 
occur, for n > m, as long as it is energetically allowed. 
It can in fact occur, from any ionic species, during the 
interaction, but it is mostly from the neutral, where all 
of the population resides at t = 0, that such processes 
are expected to play an important role. In our equations 
above, those processes are represented by the terms in 
the right hand side which contain the population Nq of 
the neutral, involving also multiphoton cross sections of 
the appropriate order. 

The values for direct n-photon m-electron ejection 
cross sections represent new territory. There is, however, 
a basis for conjecturing reasonable values. We argue that 
the cross section for n-photon, m-electron ejection (with 
m < n) in a given atom, is of the same order of magnitude 
as an n-photon one-electron ejection, for the same pho- 
ton energy [l6j . The reasoning here rests on the follow- 
ing properties of a generalized cross section: (a) The n- 
photon transition matrix element from the ground state 
connecting to n electrons in the continuum requires no 
correlation, (b) As a consequence, the summation over 
intermediate states is dominated, in both cases, by single 
electron matrix elements, (c) All intermediate states are 
in the continuum, with no intermediate resonances of any 
significance. E.g., for 2-photon, 2-electron ejection from 
an initial state Is 2 , one would have 

E (k'p, kp\ f i + f 2 \kp: Is) (kp; lsj f i + f 2 | Is 2 ) 
£ kp - E ls 2 -huj ' U 

with f.j the electronic coordinates, and the summation 
extending over all allowed 1-electron excited states. This 
is non-zero even for non-interacting electrons, reducing to 
the product (k'p\ ?2 | Is) (kp fi Is) of 1-clcctron matrix 
elements. One example of applicability of this conjecture 
is found in the cross section of two-photon double ioniza- 
tion of Helium which turns out to be roughly equal to the 
two-photon one electron ejection [l7| : while in contrast 
the cross section for one-photon two-electron ejection is 
smaller than the one-photon one-electron ejection by a 
factor of about 50. The reason of course has to do with 
the fact that ejecting m electrons by n-photon absorp- 
tion, for n < m, is impossible without correlation. A sec- 
ond example is 4-photon double ionization of Carbon Q ■ 
Adopting the above conjecture, we obtain single-photon 
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and two-photon cross sections through a calculation and 
higher order cross sections through a procedure of scaling 

am 

The result of a calculation including these direct higher 
non-linearity channels is shown by solid lines in Fig. [T] 
Their presence does not alter significantly the ion yields 
for the 30 fs pulse, but for the 5 fs pulse it does in- 
crease rather dramatically the yields of the higher charge 
species, beginning with Ne +4 . There is a clear physical 
interpretation of these results. The effect of the direct 
channels is more pronounced for the shorter pulse, be- 
cause the sequential channels do not have as much a 
chance to drain the neutral — from which the direct 
channels originate — as they have in the rising wing of 
the longer pulse. And the direct channels, being of higher 
non-linearity, will dominate only if exposed to higher in- 
tensity, provided there is still enough population left in 
the neutral; which is the case for the shorter pulse. For 
the same reason, the yields of lower charge ions (up to 
Ne +4 ) are not affected much because, in their case, even 
the direct channels are of low order. 

A first conclusion at this point is that, if the temporal 
shape of the pulse were deterministic, shortening its du- 
ration would enhance the contribution of the direct chan- 
nels. It is, however, well known that, at least for the time 
being, the pulses of the FEL sources are not deterministic 
but exhibit strong intensity fluctuations. Qualitatively 
speaking, they can be considered chaotic @. Specifically, 
the overall pulse envelop is known to contain spikes of 
random height and duration, which can be as short as 5fs. 
This means that the atom is indeed exposed to spikes of 
duration sufficiently short to enhance the direct channels. 
The implications for the theory is that Eqs. ([l) become 
stochastic, owing to the stochastic nature of the inten- 
sity. The simplest way to account for this, is to recall 
that n-photon ionization, within LOPT, is proportional 
to the 7i-th order intensity correlation function [l8| , 
which for a chaotic field is given by n\F n , where F(t) is 
the average intensity (flux) . We could thus replace F n by 
n\F n in Eqs.(l), which amounts to effectively increasing 
the n-photon ionization cross section by a factor of n\. 
Obviously, this favors processes of higher order. 

This procedure would be rigorous if we had a single 
n-photon process and in addition the field were truly 
chaotic. What we have, however, is a set of differential 
equations coupling processes of various orders. Multi- 
plying the cross sections by n\ and using a deterministic 
pulse amounts to a decorrclation approximation, valid if 
the ionic populations do not change appreciably on the 
scale of the intensity fluctuations. An alternative but rig- 
orous approach is to introduce an appropriate stochastic 
model for the radiation, and solve the differential equa- 
tions for a sufficiently large number of realizations of the 
radiation, taking in the end the average over such runs. 
This corresponds exactly to the manner experimental 
data are obtained. By comparing the results with those 
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FIG. 2: (Color online) Ionization of Ne at 93 eV under chaotic 
pulses of duration 30 fs and coherence time 6 fs. (a) Ion yields 
obtained by solving Eqs. JTJ for randomly chosen F(t), in the 
presence of sequential (dashed) and sequential+direct (solid) 
channels. The presented yields are averaged over 10 4 realiza- 
tions, (b) As in (a) with volume expansion effects included. 



obtained through the above mentioned decorrelation, we 
can also assess the limits of validity of the latter. Space 
does not allow a detailed comparison of the two models, 
nor a detailed description of our stochastic modeling. 

In the following we focus on results obtained within 
the ab initio stochastic modeling of the field, relying on 
theoretical as well as experimental information about the 
FEL pulses (e.g., see Comparing the averaged 

ion yields of Fig. Uta), to those of Fig. [IJa), we note 
a dramatic increase of the higher charge species in the 
presence of chaotic light, when the direct channels are 
included. We have thus shown that intensity fluctuations 
will also enhance the direct channels over the sequential, 
even for the longer pulse of 30fs. 



To compare with experimental data on Neon |14| . we 
have performed an integration over the spatial distribu- 
tion of the radiation in an interaction volume approxi- 
mating the experimental arrangement. The result of a 
calculation for a 30 fs pulse, corresponding to that of 
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FIG. 3: (Color online) Ionization of Ne under chaotic pulses of 
duration 30 fs and coherence time 6 fs. (a) Data obtained by 
(b-d) Ion yields at three different intensities, as obtained 
by Fig. [21b). The thin (black) bars correspond to sequential 
ionization whereas the broad (colored) bars refer to the case 
of open direct and sequential ionization channels. 



the experiment, including fluctuations, is shown in Fig. 
rjljb) . As expected, the populations of the ionic species do 
not decrease beyond the saturation intensity, but exhibit 
a slow increase, due to the contribution of more atoms 
from the periphery of the interaction volume, where the 
intensity is lower than in the center. The crucial point, 
however, is that the contribution of the direct channels 
is found to dominate, for the higher species, even upon 
spatial integration. It is therefore evident that the ef- 
fect of the multiclcctron direct channels is quite robust, 
as it survives practically intact, beyond the deterministic 
pulse. Since only peaks of TOF (Time of Flight) re- 
sults, at a single laser intensity, have been given in [lij ]. 
we have estimated the relative magnitude of those peak 
heights and have plotted as a histogram the correspond- 
ing values in Fig. G2[a). The histograms of Figs. Elb-d), 
correspond to the ion yields at three different intensities, 
at and around the nominal experimental one, as obtained 
from our Fig. [2jb) . The comparison with the experimen- 
tal data demonstrates an overall reasonable agreement, 
especially for intensities ~3x 10 15 W/cm 2 , but does also 
point to the sensitivity of the results to the value of the 
intensity. In the absence of experimental laser intensity 
dependences of the ionic species, we cannot offer at this 
time a more detailed evaluation of the agreement between 
theory and experiment. 

We can nevertheless attempt to address the question 
of whether there are traces of the direct channels in the 
limited experimental data at our disposal. To this end, 
in Figs. [3jb-d) we present the relative heights due to 
the sequential channels alone (narrow black bars) as well 
as the ones when both direct and sequential channels 
are present (broad colored bars). Depending on the ex- 
act intensity, the effect of the direct channels, although 



rather small, is noticeable particularly for higher ions in 
Figs. [3]Jb,c). The reason for the small contribution of 
the direct channels is that the data have been obtained 
at or beyond the saturation intensity; not the optimal 
conditions for detecting the contribution of the direct, 
as clearly demonstrated in Figs. HJa) and O Data over 
a range of intensities, below saturation, would provide 
a much clearer picture. Finally, somewhat surprising at 
first glance, for all three intensities the yield for the first 
ion (Ne + ) in the case of sequential ionization only, seems 
to be larger than the corresponding yield in the case of 
scqucntial+direct ionization. This, however, is consistent 
with the underlying physics, since the absence of direct 
channels leads to a much slower depletion of the neutral, 
from which direct channels originate, in favor of Ne + . 
For all other ions, the situation is reversed for all three 
intensities. 

In summary, we have introduced the idea of direct 
multiple ionization, have explored its contribution in a 
specific realistic context, having incorporated all relevant 
aspects, such as intensity fluctuations and spatial radi- 
ation distribution. Reasonable agreement with existing 
but limited experimental data has been found, highlight- 
ing at the same time the need for more detailed data. The 
results in this work provide a benchmark case and guide 
for the planning of future experiments with the new FEL 
sources which are undergoing an explosive expansion. 
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